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I INTRODUCTION 

Small-capacitance superconducting tunnel junctions provide an ideal system for studying the 
interplay between the Josephson phase and the charge on the junction electrode, which are 
quantum mechanically conjugate to each other. This interplay can be probed through the 
superconductor-insulator (SI) transition IT], which is a quantum phase transition occurring at 
r = 0. The SI transition has been extensively studied in two-dimensional (2D) systems. Ex- 
periments have been carried out on granular ||2l|3l@l|3 and homogeneous films f6{U\. Theo- 
retical studies have modeled the films as 2D arrays of small-capacitance Josephson junctions 
(JJs) ll8ll9l fT0lfTnfT2t . and experiments with such arrays have also been reported II131I14I . In ID, 
however, experimental study of the SI transition in JJ arrays has been less extensive, while the 
theoretical investigation has been done dlSEllISlElIIHlIIllEOl. Experimental data on long 
and narrow films are available ||22l|2Tl. In contrast to films, JJ arrays can be fabricated with a 
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high degree of uniformity, and the parameters of interest in the theory can be measured. Further- 
more, one can design a JJ array in such a way that one of the important parameter in the theory, 
the Josephson coupling energy Ej between adjacent islands can be tuned in situ II231I24I . The 
other important parameters are the charging energy associated with the junction capacitance, 
Ec = e^/2C, and the stray capacitance of each island to the ground, Eq^ = e^/2CQ. Depending 
on the values of these parameters either superconducting or insulating behavior is expected for 
an array with infinite length. 

In long arrays, it is possible to observe a well developed Coulomb blockade ll25ll for Cooper 
pairs in the current-voltage (I-V) characteristics, even when the Josephson energy is dominant, 
Ej > Ec- Such a Coulomb blockade is extremely interesting because our usual notions about 
phase coherence in the sense of the Josephson effect do not apply. The phase of the supercon- 
ducting state is washed out by strong quantum fluctuations, and the number of Cooper pairs on 
the island becomes well defined. Nevertheless, the large Josephson coupling causes the poten- 
tial associated with one excess Cooper pair to spread out, and in this sense the single excess 
Cooper pair becomes delocalized. In the Coulomb-blockade state, the single excess Cooper 
pair can be described as a charge soliton II261I27L which is dual to the Josephson fluxon of ID 
parallel arrays. The charge-soliton length, — 2e/2nCQVc, gives the length scale over which 
the potential is screened. Here Vc is the critical voltage for Cooper-pair tunneling in a single 
junction, which is a function of the ratio Ej/Eq [(28l. For Ej > Ec, Vc is reduced exponentially 
to zero as Ej increases. This weakening of the Coulomb blockade causes Xs and we expect 
the insulating state of the array to eventually give way to superconductivity when Ej ^ Ec- In 
SecHni we describe experimental data which display this evolution of the insulating state as Ej 
is tuned in situ. We interpret the data qualitatively in terms of a theoretical model for a T = 
quantum phase transition lH). 

We have also used the ID JJ arrays to bias a single Josephson junction in order to control 
the electromagnetic environment for the single junction ll29l . In single junctions, experimental 
observation of Coulomb blockade has been considered to be extremely difficult because a high- 
impedance environment is necessary, and special care should be taken with the measurement 
leads ll25l . For this reason, thin-film resistors 11301 and tunnel-junction arrays 113111321 were 
employed for the leads, and an increase of differential resistance around V = was reported. 
In contrast to the earlier works Il301l31ll32l . our leads are tunable, and we can therefore study 
the same single junction in different environments. We show that the I-V curve of the single 
junction is indeed sensitive to the state of the environment. Furthermore, we can induce a 
transition to a Coulomb blockade in the single junction when the zero-bias resistance of the JJ 
arrays is much higher than the quantum resistance Rk = h/e^ ^26 kQ. In addition to Coulomb 
blockade, we have clearly observed a region of negative differential resistance in the I-V curve. 
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FIG. 1: One-dimensional array of small-capacitance dc SQUIDs. (a) Scanning electron micrograph, (b) 
Schematic diagram. 

The negative differential resistance appears as a result of coherent tunneling of single Cooper 
pairs according to the theory of current-biased single Josephson junctions II251I33I . Based on 
the theory, we have calculated the I-V curves numerically. The measured I-V is consistent with 
the numerical calculation. 

II EXPERIMENT 



A Sample fabrication and characterization 

The ID JJ arrays were fabricated on a Si02 substrate with electron-beam lithography and a 
double-angle-evaporation technique ll34l . The arrays are made of Al with an AI2O3 tunnel bar- 
rier, and each of the Al electrodes in the array is connected to its neighbors by two junctions in 
parallel, thus forming a superconducting quantum interference device (SQUID) between near- 
est neighbors. Figure [T] shows a scanning electron micrograph of a section of an array and the 
schematic diagram. The advantage of the SQUID geometry is that we can control the effective 
Ejhy applying an external magnetic field, B, perpendicular to the substrate. 
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where Aioop is the effective area of the SQUID loop and 4>o = h/2e = 2 x IQ-^^ Wb is the 
superconducting flux quantum. Because the geometrical inductance of the SQUID loop, Lq <^ 
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4>o/27r/to, an external magnetic field creates a phase shift so that the critical current between 
nearest neighbors is modulated in a periodic way. Here, 



is the Ambegaokar-BaratofF critical current ll35l . which is calculated from the superconducting 
energy gap, Aq (= 0.2 meV for Al), and normal-state tunnel resistance of the junction, 
Henceforth, we will refer to the lumped SQUID as an effective junction with a tunable Ej, and 
a fixed charging energy Ec = e^/2C, where C is the sum (parallel combination) of two junction 
capacitances. 

There are a couple of ways to obtain The resistance of the array divided by the number 
of junctions, A^, measured above the superconducting transition temperature, Tc (= 1.2 K for 
Al), or in a magnetic field strong enough to completely suppress the superconductivity (> 0. 1 T 
for Al), is R„ by definition. It is also possible to find 7?„ by taking the slope of the I-V curve 
at high bias, V > N{2Ao/e). The capacitance C = c^A is estimated from the junction area A, 
where the specific capacitance is on the order of 10^ fF//jm^ as we will see later in Sec. lIVIE)! 
Another important parameter of the array is the capacitance of each electrode to the ground, 
Co ~ 10 aF [i36il . which depends on N logarithmically. In Sec. [nil we will discuss three arrays 
with nominally identical junction parameters [Rn = 4.9 kQ, A = (0.4 x 0.1 i-irn^) x 2, andAioop = 
0.7 X 0.2 ium\ but having a different A^: 255, 127, and 63. 

The samples for Sec. IIVICI and iDl are single Josephson junctions biased with the SQUID 
arrays. (See Fig. |2|) The single junction, which is in the center of Fig. |2^, has an area of 
0.1 X 0.1 /jm^. On each side of the single junction there are two leads enabling four-point mea- 
surements of the single junction. A part of each lead close to the single junction consists of 
the SQUID array with A = (0.3 x 0.1 ^m^) x 2, and Aioop = 0.7 x 0.2 ijm^. The samples are 
characterized by Rn, N, and the normal-state tunnel resistance of the single junction, r„. 

B Low-temperature measurements 

The I-V curves and the zero-temperature resistance were measured in a ^He-'^He dilution re- 
frigerator at 0.02 — IK. The temperature was determined by measuring the resistance of a 
ruthenium-oxide thermometer [|37J fixed at the mixing chamber. Special care was taken to fil- 
ter the sample from high-frequency electromagnetic radiation B4il . The preamplifier stage of 
our measurement scheme was specially designed for the high resistances associated with the 
Coulomb blockade. 
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FIG. 2: Single Josephson junction biased with arrays of small-capacitance dc SQUIDs. (a) Scanning 
electron micrograph, (b) Schematic diagram. 

The temperature dependence of the zero-bias resistance (Fig. ^ was determined with lock- 
in technique at 13 Hz ll23l . The power dissipation was kept at 10^^^ W, which was just large 
enough to yield a detectable signal, and at the same time, small enough to probe the "linear" 
response. 

The I-V curve of the single junction (Fig. |9l) were measured in a four-point configuration, 
where the potential difference was measured through one pair of SQUID-array leads with a 
high-input-impedance instrumentation amplifier, and through the other pair of SQUID-array 
leads, the bias was applied and the current was measured with a current preamplifier [12911 . 
When the voltage drop at the SQUID arrays was much larger than that at the single junction, 
the single junction was practically current biased. The SQUID arrays could be measured in a 
two-point configuration (same current and voltage leads) on the same side of the single junction. 
Note that the two arrays are connected in series and that current does not flow through the single 
junction. 

Ill SUPERCONDUCTOR-INSULATOR TRANSITION 
IN ONE-DIMENSIONAL ARRAYS 

A Current- voltage characteristics and the zero-bias resistance 

Figure [3^ shows the I-V curve of the three arrays at zero magnetic field. The arrays were made 
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Magnetic field (mT) 



FIG. 3: Dependence of the I-V curves on the array length N, atT = 0.05 K. (a) The I-V curves at 6 = 
showing Josephson-like behavior and the critical current 7^. (b) The I-V curves at B = 7.1 mT showing 
the Coulomb blockade of Cooper-pair tunneling and the threshold voltage Vf. (c) The magnetic field 
dependence of 4. (d) The magnetic field dependence of V;. 
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on the same chip, in the same vacuum cycle, using masks written to the same dimensions. Thus, 
all junctions in each array should be identical. The arrays are not truly superconducting, and 
there is actually a slope on the "zero-voltage" branch of the I-V curve, which gives a finite 
resistance. Furthermore, the observed "critical currents," i.e., the first local current maximum at 
^ 0.03 mV, are only 1% of the classical Ambegaokar-Baratoff value. This critical current shows 
a clear dependence on the array length. The longer the array, the lager the critical current, 
indicating that superconducting behavior is favored in the longer array. As Ej is suppressed 
below EjQ with an externally applied magnetic field, the measured critical current of each array 
is reduced, and the resistance on the "zero-voltage" branch increases. Figure |3l: shows the 
magnetic-field dependence of the critical current. In the neighborhood of 5^ = 5.8 mT, the 
curves in Fig. |3l: cross one another, so that for B > Be, the longer the array, the smaller the 
critical current. 

Figure |3b shows the I-V curve of the three arrays at 5 = 7.1 mT {>Bc). Here we see a 
new type of behavior which is dual to the B < Be behavior. The I-V curve is characterized by 
a zero-current state for voltage below a threshold voltage, where the array switches to a finite 
current state. The magnetic-field dependence of the threshold voltage for B > Be is shown in 
Fig. |3jl for the three arrays. We see that the longer the array, the larger the threshold voltage, 
indicating that insulating behavior is favored in the longer array. 

Figure m shows the temperature dependence of the zero-bias resistance, Ro(T), taken at 
the same magnetic fields (same Ej) for two arrays of different length: = 255 (solid) and 
63 (dashed). Each set of the curves shows qualitatively similar behavior. At zero magnetic 
field, as the temperature is lowered Rq decreases to a value which is temperature independent. 
As the magnetic field is increased, the resistance of this "flat tail" increases, until it reaches a 
critical value, where Ro{T) curves make a sharp turn to increasing resistance as T ^ 0. Further 
increasing of the magnetic field drives the array into the insulating state, where Rq increases 
rapidly as T ^ 0. 

If we examine the bottom two curves in each set of Fig.|4l we can see that at high tempera- 
tures, the 63-junction array has a smaller resistance than the 255-junction array, as expected for 
a classical resistor. However, at low temperatures, the resistance of the 63-junction array be- 
comes larger than the 255-junction array. This increasing of the resistance for shorter arrays is 
a clear sign that quantum fluctuations are responsible for the measured resistance [IJ. The open 
circles in Fig.|4|indicate the crossing points where Rq is the same for two different lengths at the 
same magnetic field. This crossing point moves towards T = as the magnetic field is tuned to 
the critical point . At the critical point, the T ^ resistance due to quantum fluctuations, is 
independent of the array length. 

We have seen a magnetic-field-tuned transition from Josephson-like behavior to Coulomb 



7 



10« 




1Q4I ^ ^ ^ ^ 1 
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FIG. 4: Zero-bias resistance vs. temperature for two arrays having different number of junctions, A'^, 
but otherwise identical parameters. The set of soUd curves are for A'^ = 255 and the dashed curves are 
for N = 63, taken at the same magnetic fields between and 7 mT. The open circles show where the 
measurements on the two arrays at the same magnetic field cross. At the magnetic field where Ko = Kq, 
the two arrays have the same T — > resistance, which is presumably independent of A^. 
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FIG. 5: Typical path or time history of a one-dimensional Josephson-junction array in (a) the supercon- 
ducting phase and in (b) the insulating phase, respectively. The orientation of the arrows indicates the 
phase angles of the superconducting order parameter on the metallic elements connected by the Joseph- 
son junctions. 

blockade, which may be called a superconductor-insulator transition, through the 1-V curve 
(Fig.© and the temperature dependence of i?o (Fig-Si- Moreover, the sharpness of the transition 
is strongly influenced by the length of the array. We can find qualitative explanation for this 
length dependence in a theoretical model of a quantum phase transition, which will be discussed 
in the following subsection. 

B Mapping to the XF model 

The SI transition can be described in an elegant theoretical framework as a quantum phase 
transition In these models, one can describe how a T = property of a macroscopic quan- 
tum system with many degrees of freedom, will change as the complementary energies in the 
Hamiltonian of the system are adjusted. One can calculate the linear response, which in our 
case is the zero bias-resistance Rq, resulting from quantum fluctuations of the degrees of free- 
dom. Within this framework, our ID quantum system of Josephson junctions is mapped to 
the classical XY model of (l-i-l)D as sketched in Figs.|5t and|5j5, the extra dimension being 
imaginary time, ih/kBT . Note that the role of temperature for the quantum system is to set 
the "size" of the system in the imaginary-time dimension. The (l-i-l)D classical XY model 
exhibits a Berzinski-Kosterlitz-Thouless phase transition B81I39L from a disordered state (free 
vortices) to an ordered state (bound vortex pairs) as the strength of the dimensionless coupling 
constant, is increased. The quantum fluctuations of the phase of the superconducting wave 
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FIG. 6: Zero-bias resistance taken at T = 0.05 K as a function of the dimensionless coupling constant 
Kq = {Ej /^EcqY^^ for three arrays having different number of junctions. 



function are thus described in terms of vortices, and in the insulating state (large quantum fluctu- 
ations) corresponds to the free- vortex state of theXF model ||40|. The mapping to the isotropic 
XY model ifTSi can be done in the limit Q 3> C which results in Kq = [Ej/SEQy^^, where 
= e^/2Co is the charging energy associated with the stray capacitance of each electrode. 
Thus in this mapping, the junction capacitance, C is neglected. 

Figure m shows a plot of Rq measured at the lowest temperature, T = 0.05 K, vs. = 
[Ej/SEcqY^^ for three arrays with A'^ = 255, 127 and 63. We see in Fig.|6lthat the three curves 
for different lengths cross at nearly the same point, Kq = 0.027. To the left of this crossing point 
we have the insulating state, where the resistance is larger for the longer arrays. To the right 
of this crossing point, we have the superconducting state, where the resistance is smaller for 
longer arrays. The arrays in the experiments have C ^ Co, and thus we can not directly apply 
the theory. In our earlier work ll23l we postulated that the effect of C ^ Co could be accounted 
for by "course graining" to the scale A = (C/Co) , which would result in a coupling constant 
J = (Ej/AEc) Figure 4 of Ref. ll23l shows that for this choice of the dimensionless coupling 
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constant, the Rq vs. J curves do not cross at the same point, but in the region 7 G {0.49,0.55}. 
Choi et al. 11201 have made a theoretical analysis of the role of a finite junction capacitance C 
by treating A as a small parameter. They found that the transition point should approach the 
limiting value KQ = 2/n = 0.64 when extrapolated to large A. The fact that the curves cross at 
one point in Fig. |6l would suggest that = {Ej/SEcqY^'^ is indeed the correct parameter for 
the transition. However, the experiment does not support the conclusions of Choi et al. in that 
the experimental critical point ^ 0.03 <^ Kq = l/n. 

As we have seen earlier in Fig. |4l the resistance at T > 0.6 K is almost proportional to 
the array length, and in this sense, the arrays behave like classical ID resistors. At lower 
temperatures, however, large deviations from this classical behavior occur. We can qualitatively 
understand these observations in the context of the ( 1 + 1)D XY model by considering the finite- 
size effect, which plays an important role in a real experiment. The zero-bias resistance of the 
array is determined by quantum fluctuations of the superconducting phase, which are described 
by the vortices in the (l-i-l)D XY model. As the temperature is lowered, the system becomes 
larger in the imaginary-time dimension, and at low enough temperatures, the system size in the 
real- space dimension, or the array length, determines the energy for free- vortex formation. The 
energy increases with increasing system size, and thus the probability of free-vortex formation 
is reduced in a longer array. This means that in a longer array, the superconducting state is 
favored. If we assume that the zero-bias resistance of the array is proportional to the probability 
of free- vortex formation in the isotropic XY model with the area of A^^ (i.e., units in real space 
and N units in imaginary time), we obtain 

Rq ~ A^2-7r/fo (4) 

by neglecting any renormalization effects El- For > KQ = 2/n, Rq increases for decreasing 
A'^, as observed in Fig. [Hon the superconducting side of the transition. 

IV COULOMB BLOCKADE IN SINGLE JUNCTIONS 
A Theory for current-biased single Josephson junctions 

The Hamiltonian of a single Josephson junction in an environment with sufficiently high impedance 
is written as 

H = —-Ejcos(^, (5) 

where Q is the charge on the junction electrode, C is the capacitance of the junction, Ej is the 
Josephson energy, and {|) is the Josephson-phase difference across the junction. The charge Q 
and h(^/2e are quantum mechanically conjugate valuables, and a set of the eigenfunctions are 
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Bloch waves of the form 

y\f{(^) = m((^) exp(%/2e), (6) 
where q is called quasicharge and u{^) is a periodic function, 

a((|) + 27r) = (7) 

The energy eigenvalue E plotted as a function of q has a band structure, and in all the allowed 
bands, it is 2e periodic. An example of the energy diagram for Ej/Ec = 0.2, where Ec = e^/2C 
is the charging energy, is shown in Fig. |7^. Under constant current bias /x, in the absence of 
quasiparticle or Cooper-pair tunneling, q increases uniformly in time according to 

so that the state of the system advances toward higher q within a given band as time goes on. 
The average voltage is given by 

(V) = En..,</)^^. (9) 

where ib is the band index and P{ih,q) is the probability that the system is in the state {ib,q)- 
The probability P{ib,q) can be calculated by solving a set of coupled differential equations of 
the form 

= l^A{ib,q4,q')P{i'b,q') =0, (10) 

where the matrix element A(/i, q, i'^, q') describes the rate of transition between the states (Z^,, q) 
and {i'b^q')- The dominant process for A(z/„^,z^,^') depends on the magnitude of h- 

An example of the theoretical I-V curve is shown in Fig.lVj). For sufficiently small 4 (re- 
gion CB), the dominant process is stochastic quasiparticle tunneling, where q changes by e. 
This tunneling always occurs along the energy parabola q^ /2C, such that it changes by or 
+ 1 if the initial state is in the lowest band {iy = 1) and by ±1 for all the other initial states 
{ib > 2) ll33l . The rate for the quasiparticle tunneling is given by 

, , AE/e^Rnr, 
^'^'- exp(Ai;/W-l - 

where Rqp is the quasiparticle resistance and AE is the difference in energy between the initial 
(ib.q) and final {i'b.q') states, 

AE = E{ilq')-E{ib,q). (12) 

At sufficiently low temperatures, the tunneling with > is extremely unfavorable, and the 
I-V curve is highly resistive (Coulomb blockade). For larger 4 (region BO in Fig. |7j5), the 
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FIG. 7: (a) Energy diagram and (b) theoretical current-voltage characteristics for a single Josephson 
junction with Ej/Ec = 0.2 and Rqp = 200 {h/n^e^) at ksT/Ec = 0.2, where Ej is the Josephson energy, 
Ec = e^/2C is the charging energy, Rqp is the quasiparticle resistance, and ksT is the thermal energy. 
(The energy diagram depends only on Ej/Ec.) 
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quasicharge is frequently driven to the boundary of the Brillouin zone, q = e, then taken to —e 
as a Cooper pair tunnels (Bloch oscillation). This process decreases (V), and as a result, the I-V 
curve has a region of negative differential resistance, or "back bending" in the low-current part. 
For still larger 4 (region ZT in Fig.|7j5), Zener tunneling becomes important, and (V) increases 
again. In Zener tunneling, no quasicharge is transferred but the state of the system jumps from 
one band to another as it passes by the narrow gap between the bands. The probability of Zener 
tunneling from band to ih + 1 or vice versa is given by 

n {AEf e 



exp 



(13) 



8 ibEc hh 

Following Ref. ll42l which takes into account the above tunneling processes (quasiparticle. 
Cooper-pair, and Zener), we have calculate the I-V curve numerically. The parameters for 
the calculation are Ej/Ec, hgT /Ec, and a = h/n^e^Rqp. The current and the voltage are in 
units of e/RqpC and e/C, respectively. 

As we have seen in Fig.|35, a typical I-V curve consists of three regions, so that it is charac- 
terized by the local voltage maximum, or blockade voltage Vb, and the local current minimum, 
or crossover current 1^. (Here, we have to mention that the back-bending feature is smeared out 
if kgT /Ec or a is increased considerably.) Analytic expression of Vb and has been obtained 
theoretically for limiting cases ll33l . The value of Vb is a function of Ej/Ec, and given by 



^ _ 0.25 e/C for£y/£c<l, 
^"^^ 5o/e ioxEj/Ec->l, 



(14) 



as r ^ 0, where 



5o 



C 



27^2 



1/4 



(I) 



3/4 



exp 



Ec 



1/2 ■ 



(15) 



is the half width of the lowest energy band. As for /„, 



RqpC 



1/2 



(16) 



is expected for a <^ [Ej/EcY ^ 1 and T ^Q, where 



Iz = 



n eE; 



(17) 



8 hEc 

is the Zener breakdown current. Note that I^r is much smaller than Iz- When we compare 
our experimental results with the theory, we need theoretical prediction for finite kgT /Ec, and 
arbitrary Ej/Eq and a. For this reason we have done the numerical calculation. The measured 
Vb and /cr will be compared with the calculation in Sec. IIVIDI 
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Voltage (2.5 |lV/div.) 

FIG. 8: Current- voltage characteristics of a single Josephson junction biased with thin-film resistors at 
several temperatures. From left to right, T = 0.06, 0.08, 0.11, 0.17, 0.23, 0.28, and 0.34 K, respectively. 
The origin of the voltage axis is displaced for each curve for clarity. 

B Earier experiments on single junctions 

In order to observe the Coulomb blockade in a single junction experimentally, the electromag- 
netic environment for the junction, or the measurement leads connected to the junction, should 
have a high impedance ll25l . For this reason, thin-film resistors (NiCr alloy, AuPd alloy, Cr, and 
layered Ge/Pd) ll30ll and tunnel-junction (AI/AI2O3/AI) arrays |l3Tl|32j| were employed for the 
leads to bias a single junction (AI/AI2O3/AI). 

FigurelUshows the I-V curves of a single Josephson junction biased with thin-film resistors 
at several temperatures. A clear Coulomb blockade is seen at T < 0. 11 K, and a "back bending" 
is also visible. In similar samples, high-frequency (/ = 0.4 — 10 GHz) irradiation induced steps 
in the dc I-V curve at / = ±2e/ (not ±e/), which can be explained as a phase locking of the 
externally applied signal to the Bloch oscillations [30|. 

Geerligs used 2D tunnel-junction arrays for the leads ll3Tl . He claimed that in the normal 
state {T > Tc), a Coulomb blockade was visible in the single-junction I-V curve. In the su- 
perconducting state (T < Tc), however, neither the arrays nor the single junction develop clear 
charging effects, and the single junction showed a classic hysteretic I-V curve with an ordinary 
supercurrent. We believe that his arrays were in the superconducting side of the SI transition, 
and did not have high enough impedance below T^. 

Shimazu et al. biased a single junctions with ID tunnel-junction arrays, and reported an 
increase of differential resistance around V = 0m the normal state ll32l . In the superconducting 
state, they measured the zero-bias resistance rather than the I-V curve. The zero-bias resistance 
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in the superconducting state was higher than the normal-state resistance for some single junc- 
tions, which suggests the existence of a Coulomb blockade even in the superconducting state. 
However, this increase of the zero-bias resistance could also be due to simple quasiparticle 
tunneling, which is independent of any Coulomb blockade effects for Cooper pairs. 

C Single Josephson junctions biased with SQUID arrays 

We have employed ID arrays of dc SQUIDs for the leads to bias a single Josephson junc- 
tion ll29l . The advantage of this SQUID configuration is that in contrast to the earlier experi- 
ments in Sec. lIVIBl the impedance can be varied in situ by applying an external magnetic field 
at low temperatures. (See Fig.lH) Thus, we can tune the electromagnetic environment for the 
single junction over a wide range. In Figs. |9l and [lOl we show some results on a sample with 
r„ = 17 ktl, Rn = 1.4 kQ., and = 65. The I-V curves for the single junction at several normal- 
ized magnetic fields, cp = 5Aioop/^o> are shown in Fig.|9l As cp is varied, the I-V curve develops 
a Coulomb blockade. We emphasize that the Josephson energy of the single junction is inde- 
pendent of cp, because it does not have a SQUID configuration and the field (p4>o/Aioop < 7 mT 
applied here is much smaller than the critical field for Al films 0.1 T). The electromagnetic 
environment for the single junction (the SQUID array), however, is strongly varied with cp. The 
behavior of the single junction demonstrated in Fig. |9l does not result from the magnetic-field 
influencing the single-junction I-V curve, but rather from an environmental effect on the single 
junction. This experiment demonstrates in a direct way that the single-junction I-V curve is 
indeed sensitive to the electromagnetic environment. 

The I-V curves of the two SQUID-array leads connected in series at cp = 0.43 (Rq = 0.61 MQ), 
0.46 (Rq = 3.2 MQ), and 0.49 (Rq = 43 MQ) are shown in Fig. [JOl The I-V curves of the 
leads are nonlinear, and in general the SQUID array cannot be described by a liner impedance 
model i24J. However, we may characterize the environment by their Rq. Coulomb blockade 
is visible only when Rq S> Rk, which is consistent with the theoretical conditions for the clear 
observation of Coulomb blockade in single junctions [il43ll . For an arbitrary linear environment 
characterized by Zg(a)), Re[Ze((o)] ^ Rf- is required for the Coulomb blockade of quasiparticle 
tunneling and Re[Ze(a))] ^ Rk/4- for that of Cooper-pair tunneling ll43l . It is interesting to 
note that at cp = 0.46 (labeled "b"), the I-V curve of the leads is still "Josephson-like" (differ- 
ential resistance is lower around V = 0), while that of the single junction is already "Coulomb- 
blockade-like". This feature becomes more distinct in samples with larger A'^ B4l . 

D Comparison with the numerical calculation 

The region of negative differential resistance seen in Fig.|9lwhen Coulomb blockade is well de- 
veloped, is related to coherent tunneling of single Cooper pairs according to the theory ll25l 
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FIG. 9: Current-voltage il-V) curves of the single junction at T = 0.02 K for a sample with r„ = 17 kH, 
Rn = 1.4 and N = 65. From top left to bottom right, the normalized magnetic field cp = BAioop/<I>o is 
increased from 0.43 to 0.49 in steps of 0.01. The origin of each curve is offset for clarity. For the labeled 
curves, the I-V characteristics of the leads at the same cp are shown in Fig.fTOl 
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FIG. 10: Current-voltage curves of the two SQUID-array leads connected in series at T = 0.02 K for the 
same sample as in Fig.|9l From a to c, cp = BAioop/<I>o is 0.43, 0.46, and 0.49, respectively. 

l33l of a current-biased single Josephson junction in an environment with sufficiently high 
impedance. Following Ref. |42|, we have calculated the blockade voltage V}, numerically as 
a function oiEjjEc ll45l . The measured Y\, (See Fig.|9l) for the samples having nominal junc- 
tion area of 0. 1 x 0. 1 ^um^ is compared with the numerical calculation in Fig.[TTl The boxes and 
circles represent the samples biased with thin-film resistors fSec. lIVIBll and with SQUID arrays 
(Sec. lIVICb . respectively. For the sample shown in Fig.|9l we used the data at cp = 0.49 (curve c) 
in order to obtain V^,. At cp = 0.49 the voltage drop at the SQUID arrays is 10^ times larger than 
that at the single junction, and the single junction is therefore considered to be current biased. 
Compare the voltage scale of Figs.|9land[l0l We calculated Ej from r„, Ej = /zAo/8e^r„. For 
Ec, we employed c., = 130 fF//jm^, and with this value the experimental data, especially those 
for the samples biased with SQUID arrays (the circles in Fig. fTTT) . agree with the numerical 
calculation. 

Actually, a smaller value, Cs = 45 ± 5 fF//jm^ ll46l , which was obtained for the junctions with 
3 X 28 jjm^ and 7 x 54 pm^, has been frequently employed ||23l|24l|30l|32l|36|. Our apparently 
large Cs may be partly explained by distributed capacitance of the SQUID arrays or by residual 
environmental effects. We also note that the uncertainty in Cs seems to be large when the 
junction area is on the order of 0.01 fjm^ or smaller. For example, Fulton and Dolan measured 
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FIG. 11: Blockade voltage Vb divided by e/C as a function of Ej/Ec- From top to bottom, the curves 
represent the numerical calculations for normalized temperatures ksT /Ec = 0, 0.02, 0.05, 0.1, 0.2, and 
0.5, respectively. 

samples with three junctions that share a common electrode, and obtained 0.20 — 0.23 fF for 
(0.03 ± 0.01 ^m)2 X 3 ED, i.e., c, = 42 - 192 fF/^um^. Geerligs et al. reported c, ^ 110 fF/^m^ 
for two-dimensional (190 x 60) junction arrays with the areas of 0.01 or 0.04 /jrn^ ifTSl . More 
recently, Penttila et al. studied resistively shunted single Josephson junctions with the area of 
0.15 X 0.15 fim^ f4B\. The estimated C of their eight samples ranged between 0.8 and 6.6 fP, or 
c, = 36 - 293 fF/pm^. 

When C is determined, it is possible to estimate Rqp of the samples from lev We plot the 
measured 1^ for the single junctions biased with SQUID arrays (Sec. IIVICI) as a function of 
Ej/Ec together with some theoretical curves based on our numerical calculation in Fig. [121 We 
obtain i?qp = 10° - 10^ Ma, or i?qp/r„ = 10^ - lOl 

V CONCLUSIONS 

One-dimensional (ID) arrays of small-capacitance SQUIDs undergo a sharp transition, from 
Josephson-like behavior to the Coulomb blockade of Cooper-pair tunneling, as the effective 
Josephson coupling between nearest neighbors is tuned with an externally applied magnetic 
field. We have shown how length scaling of the zero-bias resistance of the array can be used 
to probe the superconductor-insulator quantum phase transition. The observed non-classical 
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FIG. 12: Crossover current /^r vs. Ej/Ec- The curves represent the numerical calculations for C = 
1.3 fF, and from top to bottom {a = h/K^e^Rqp,kBT /Ec) = (IQ-^ ,0), (10-2,0.3), (10-2,0.5), (10-^,0), 
(10-^0.3), (10-3,0.5), (10-^,0), (10-4,0.3), and (10-^,0.5), respectively. 

dependence of the zero-bias resistance on the length of the array, where the zero-bias resistance 
decreases with increasing length, can be supported qualitatively with a theoretical model which 
maps a ID quantum system to the (1 + 1)D classical XY model. 

We have also used the SQUID arrays as a tunable electromagnetic environment for a single 
small-capacitance Josephson junction, and demonstrated how the Coulomb blockade of Cooper- 
pair tunneling is induced in the single junction. When the Coulomb blockade is well developed, 
the measured current-voltage curve is consistent with the numerical calculation for a current- 
biased single Josephson junction. 
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